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Abstract 

We study the onset of synchronization in lattices of hmit cycle oscillators with 
long-range coupling by means of numerical simulations. In this regime the critical 
coupling strength depends on the system size and interaction range reflecting the 
non extensive behavior of the system, but an adequate scaling removes the depen- 
dency and collapses the long-range synchronization curves with the one resulting 
from a system with uniform coupling. Two descriptions are considered, the Ku- 
ramoto or phase description model and a Cell Dynamical System model for phases 
and amplitudes. Oscillator death is observed in the second approach. 



1 Introduction 

Coupled limit-cycle oscillators have been extensively used to model the col- 
lective behavior of systems in physics [1], chemistry and biological sciences 
[2,3]. The existence of a self organized periodic motion is a direct consequence 
of the energy exchange between the system and its surroundings keeping the 
system far from equilibrium. The mathematical description of those systems 
relies on differential equations, or its discrete versions, since no Hamiltonian 
exists, and the appearance of a macroscopically synchronized state has been 
treated as non-equilibrium phase transition [3,4]. The coupling among oscil- 
latory units is usually considered to be uniform and of infinite range or with 
first neighbors only. The reason for these choices is simplicity, uniform cou- 
pling makes possible analytical calculations, and nearest neighbor coupling is 
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a computationally efficient form of short range coupling. Here we would like 
to consider a long range distance dependent coupling of the form r~", where 
r is the distance between oscillators and a > defines the range: a = cor- 
responds to uniform coupling, and a — > oo to nearest neighbor coupling. For 
this we consider two different approaches: the Kuramoto model for the phases 
and a cell dynamical system (CDS) model for phases and amplitudes along 
the limit cycle. We observe that, as in equilibrium systems, the critical value 
for synchronization will depend on the system size and on a. For both models 
we show that a scaling proposed for equilibrium nonextensive spin systems 
works well for the synchronization process. 



2 Long Range Coupling in Equilibrium Spin Systems 



It well known that, depending on the dimensionality d of the system, inter- 
acting potentials of the type J/r~" may lead to nonextensive energies. This 
problem is usually avoided in the case of uniform or mean field coupling, for 
which a = 0, by the redefinition of the coupling term with the division by 
the size of the system, N. This correction may be extended to other values of 
a in the nonextensive regime with the definition of a N*{a,d) that may be 
calculated from the requirement that the total energy of the system must be 
finite. The equilibrium properties of chains of ferromagnetically coupled spins 
have been considered in the nonextensive regime for the Ising model [5,6] and 
XY model [7,8], the finiteness condition for the energy leads to a common 
correction factor 

/ A/1-''/'' - 1 

N*{a,N)= / dr/-V-" = — . (1) 

^ 1 — a/a 



The extensive behavior is recovered with the definition of the scaled interac- 
tion strength J' = J/{2°'N*). For a > d, (1) provides a correction for finite 
size effects, and ior a < d it removes the dependency of the critical coupling 
strength on a and N. In the limit N ^ oo, N* — N for mean-field cou- 
pling (a = 0) and, A^* = InA^ for a — d. The behavior of thermodynamical 
quantities like magnetization and energy as a function of temperature T for 
different values of N and a are found to collapse with the mean field curve 
when plotted as a function of T/{N*2°'), indicating that this regime contains 
all the information regarding the nonextensive behavior. 
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3 Kuramoto Model 



One of the simplest models for systems of coupled limit-cycle oscillators is 
the Kuramoto model for the time evolution of the phases. As shown in [4], it 
possible to reduce any system of weakly coupled limit cycle oscillators near 
the Hopf bifurcation, to the set of equations 

N 

e^ = uj, + J2^M-9i), {i = l,2,...,N) (2) 
j 



where 6i and Ui are the phase and natural frequency of the ith oscillator 
respectively. The coupling term Tij{6i — 6j) is a periodic function that may 
be obtained from the original equations of motion, regardless of the dimen- 
sionality of the oscillating field. The case of uniform, or mean field, coupling 
has been considered with detail by Kuramoto [4] for the simplest form of the 
coupling term, 

Tij ^ N-^Ksm{ej -Oi). {K>Q) (3) 



Kuramoto considered the possibihty of a collective oscillatory motion when 
the natural frequencies are picked from a given distribution g{u!), as would 
happen in a real system. As the coupling modifies the frequencies, acting as 
an external oscillatory force, it is necessary to define the asymptotic frequency 
of the ith oscillator after the transient regime as 

(1;,= hm ^[e,(t + T)-e,{t)]- (4) 

1 — >oo _Z 



The existence of a collective behavior with phase and frequency entrainment 
depends on the couphng strength and on the choice of g{uj) [4,9]. In what 
follows we consider only the case where g{ui) is Gaussian. Kuramoto charac- 
terized the degree of synchronization in the asymptotic regime by means of two 
order parameter like quantities: R for the frequencies, and a for the phases. 
Frequency entrainment may be quantified simply by analyzing the formation 
of frequency clusters and defining 

TV 

lim -f, (5) 



where A^^ is the size of the largest cluster of oscillators with a common asymp- 
totic frequency a). The complex quantity 



1 ^ 

(7= — ^exp(i^j) (6) 
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is appropriate for the phase ordering as it measures the concentration of phases 
on a certain value. R and a behave as equihbrium order parameters in the sense 
that, for K below some critical value K^, oscillations are incoherent leading 
to R and a—0. For strong enough coupling, or K > Kc, collective oscillation 
takes place and i?, cr > 0. 

For the case of uniform coupling, phase and frequency entrainment appear 
simultaneously at the critical value 

Kc = — (7) 

where ujq is the average value of the distribution of natural frequencies [4] . 
The long range version of the coupling term is written as 

r,,- = ^ sin(^,- -e,), K> = 1, 2, 3 . . .), (8) 



and was considered in the case of oscillator chains in [10]. In Hamiltonian 
systems the limit of non-extensivity can be found by requiring the upper value 
of the energy to be finite as ^ oo. Limit cycle oscillators are open systems 
and energy is a meaningless quantity so, instead, we establish the finiteness 
condition for {l/N)J2i0i in the asymptotic regime. This condition may be 
established by requiring that the upper value of the coupling term divided by 
N must be finite, as in the mean field coupling case. With this, we end up with 
a condition identical to (1), but, as will be clear from the graphs below, the 
factor 2" does not seem to be present in this case. For any value of d and a, 
it is then possible to incorporate the dependence on a and N of the coupling 
term into the value of K with the definition of a scaled coupling strength 

K*{a,N) = KN*. (9) 



In order to numerically study the non-extensive region, we used two dimen- 
sional square lattices of coupled oscillators with natural frequencies sorted 
from the Gaussian distribution g{u;) — (l/-\/27r) exp(— a;^/2) and periodic 
boundary conditions. Equation (2) with the coupling (8) was then integrated 
using the Eulcr method with dt = 0.1. The first 3000 iterations where dis- 
carded as a transient, and after 6000 iterations past the transient, we cal- 
culated the asymptotic frequencies as defined in (4). We used the Hoshen- 
Kopelman method [11] to analyze the frequency pattern, two sites were con- 
sidered to belong to same cluster if the frequency difference between them was 
less than 27r/T, where T is the time elapsed after the transient. R was then 
associated to the largest cluster, and averaged over 10 initial conditions for 
the natural frequencies. 
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The non-extcnsivity is evident if we consider the behavior of {R) as a function 
of K for a < d. For fixed a, the critical value of K^. decreases with increasing 
N as can be seen in Fig. (l.a), as more terms contribute to the coupling. The 
shift in Kc may be corrected with the use of (1). In fact, if we plot {R) as a 
function of K*, there is no dependence on the system size, the curves collapse 
and a unique value of K* appears. Figure (l.b) shows the behavior of (R) as 
a function of K* for fixed a. The dependence of Kc on a is evident if we 




Fig. 1. Frequency entrainment in square lattices with N Kuramoto oscillators inter- 
acting through a long range coupling K/r^'"^^. (a) The non-extensivity appears in 
the N dependency of Kc, the critical value for synchronization, (b) The corrected 
coupling intensity, K* , absorbs the non-extensive part of the coupling term leading 
to a synchronization curves independent of the system size. 




Fig. 2. Frequency entrainment in square lattices with 1024 Kuramoto oscillators 

interacting through a long range coupling K/r^. (a) The non-extensivity appears 
in the a dependency of Kc, the critical value for synchronization, (b) The corrected 
coupling intensity, K* , absorbs the non-extensive part of the coupling term leading 
to synchronization curves independent of a The curve for a = correspond to 
uniform or mean-field coupling. 

plot (R) for fixed N and different values of a. Figure (2. a) shows the curves 
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for a = 0.25, 0.50, 0.75, and (mean-field coupling without the correction 
factor 1/A^). It is important to notice that we did not use the factor 2" as 
predicted for the spin systems. Consistently with dependence on system size, 
as a decreases, more terms are included in the coupling, leading to smaller 
values of Kc- Again, (1) fixes the dependence on a providing a perfect collapse 
even with the mean-field curve as can be seen in Figure (2.b). Notice that the 
value of K* in Figures (1) and (2) is in good agreement with (7), our frequency 
distribution leads to = 1.59. 



4 Cell Dynamical System Model 

A more complete description of the system includes information about the 
amplitudes by means of the definition of an oscillatory field 

z{n,t) — p{n,t) exp[i4>{n,t)]. 

This is an important aspect since oscillator death is possible as we effectively 
increase the coupling strength with the variation of system size and coupling 
range [12]. 

The spirit of cell dynamical modeling is the discretization of time and space 
from the start, as opposed to formulating a continuous model and then further 
adopting some discretization scheme. In this sense, (n, t) defines the position 
of the n-th cell at time t, p and 4> are then the amphtude and the phase of 
the field at given cell. The time evolution of the system is then given by the 
repeated application of the map 

z{t, t + l)= F[z{t, t)] + Dij{F[z{j, t)] - F[z{t, t)]} (10) 

where the local dynamics is regulated by F\z(t): A] = A exp(ici;) , ^^^l 

^ J L V J y\ ^^i+|^(t)|2(^2_i) 

This form of F corresponds to an asymptotically stable limit cycle of unit ra- 
dius for all A > 1 [13]. Here the parameter A is a measure of the attraction to 
the limit cycle. The long range coupling is implemented through the definition 
of Dij^Do/rf^. 

For observation of synchronization, we considered chains oi N — 100 oscillators 
with periodic boundary conditions, A — 2.0, Gaussian distribution of natural 

frequencies with unit mean value and variance 0.1. We iterated (10) 4000 
times, the first 1000 being discarded as a transient, and the order parameter 
a was then calculated. This procedure was repeated for different chain sizes 
and values of the decay exponent a for several values of the coupling strength 



6 



Do. Figure 3(a) shows the synchronization curves averaged over 20 samples. 
As expected, the critical value of Dq is dependent on and a. Although the 
calculation of the upper bound of the coupling term is rather different from 
the one in the Kuramoto model or spin systems, the same form of scaling is 
observed here. Figure 3(b) shows the synchronization curves of cr as a function 
of D* = Do2"A^*. 




Fig. 3. Phase entrainment in a CDS chain with oscillators interacting through a 
long range coupling Dq/v'^. (a) Smaller values of a lead smaller critical values of 
the critical coupling constant, (b) Again, the definition of the effective coupling 
constant, taking into account the nonextensive part of the coupling term provides 
the collapse of synchronization curves for different exponents. 

It is important to verify the effect of nonextensivity on the stability of the 
oscillatory state. It is well known that, depending on the amplitudes and 
frequency spread, an increase in coupling strength may lead to stabilization of 
the rest state. This has been observed for systems with mean field coupling [12] 
and for pairs of oscillators [13,14]. In order to observe this effect we increased 
the frequency spread to 0.7, then, for a fixed value of A^ and a, we varied 
Dq and A. After the 4000 iterations we verified whether the amplitude was 
non zero. The same procedure was repeated for a — 0.1, 0.5 and 0.9. Figure 
4(a) shows the A ^ Dq phase diagram indicating regions of synchronization, 
incoherent motion and death for the different values of a. The scaled phase 
diagram can be seen in Fig. 4(b). 



5 Conclusion 

The above results bring the process of synchronization even closer to an order 
disorder transition, as we directly use the scaling proposed for nonextensive 
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Fig. 4. Phase diagram showing regions of synchronization, incoherent motion and 
oscillator death for a CDS chain with oscillators interacting through a long range 
coupling Dq/t"'. Relative to each curve, right side region corresponds synchronized 
(S) motion, on the left side region oscillations are incoherent (I), and inside the 
lambda oscillators are dead (D). (a) Smaller values of a corresponds to a shift to 
lower values of Dq in the region of death, (b) The same phase diagram, now with 
the scaled coupling constant, showing the collapse of the diagrams for different 
exponents. 



equilibrium systems. Also, we see that the important quantity for synchro- 
nization and stability of the oscillatory phase is the scaled coupling constant. 
With this in mind, one can easily tune the interaction strength to obtain the 
desired regime, in the case of distant dependent coupling. Some points still re- 
main unclear. First, there is no a priori reason for the scahng of the CDS model 
to be the same as that of Hamiltonian spin systems. Second we cannot justify 
the absence of the factor 2" in the scaling of the Kuramoto model, which has 
a condition for extensivity identical to Hamiltonian spin systems. Actually we 
found by tentative and error that for the Kuramoto oscillators a better form 
of scaling should have a factor 2"/"^. Since the results for spin Hamiltonians 
and for the CDS model are for one dimensional systems, this correction would 
not affect their scahng, although we are not able to analytically justify this 
form. 



acknowledgments 



The authors thank Francisco Tamarit and Celia Anteneodo for interesting 
discussions and the Brazilian agencies CNPq, FAPERJ and FUJB-UFRJ for 
financial support. 



8 



References 



[I] D. S. Fisher, Sliding charge-density waves as a dynamic critical phenomenon, 
Phys. Rev. B 31 (1985) 1396. 

[2] L. Glass, M. Mackey, Prom clock to chaos.The rhythms of life, Princeton 
University Press, Princeton, 1988. 

[3] A. T. Winfree, The Geometry of Biological Time, Springer- Verlag, New York, 
1980. 

[4] Y. Kuramoto, Chemical oscillations waves and turbulence. Springer- Verlag, 
Berlin, 1984. 

[5] F. D. Nobre, C. Tsallis, Infinite-range ferromagnet:thermodynamic limit within 
generalized statistical mechanics, Physica A 213 (1995) 337. 

[6] S. A. Cannas, F. A. Tamarit, Long-range interactions and nonextensivity in 
ferromagnetic spin models, Phys. Rev. B 54 (1996) R12662. 

[7] F. Tamarit, C. Anteneodo, Rotators with long-range interactions connection 
with the mean-filed approximation, Phys. Rev. Lett. 84 (2000) 208. 

[8] A. Campa, A. Giansanti, D. Moroni, Canonical solution of a system of long- 
range interacting rotators on a lattice, Phys. Rev. E 62 (2000) 303. 

[9] H. Daido, Lower critical dimension for populations of oscillators with randomly 
distributed frequencies; a renormalization-group analysis, Phys. Rev. Lett. 61 
(1988) 231. 

[10] J. L. Rogers, L. T. Wille, Phase transitions in nonlinear oscillator chains, Phys. 
Rev. E 54 (1996) R2193. 

[II] J. Hoshen, R. Kopelman, Percolation and cluster distribution. I. Cluster 
multiple labeling technique and critical concentration algorithm, Phys. Rev. 
B 14 (1976) 3438. 

[12] G. B. Ermentrout, Oscillator death in populations of "all to all" coupled 
nonlinear oscillators, Physica D 41 (1990) 219. 

[13] M. Bahiana, M. S. O. Massunaga, Cell-dynamics modeling of oscillator systems, 
Phys. Rev. E 52 (1995) 321. 

[14] K. Bar-Eli, On the stability of coupled chemical oscillators, Physica D 14 (1985) 
242. 



9 



